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ABSTRACT 

We derive masses and radii for both components in the single-hned echpsing 
binary HAT-TR-205-013, which consists of a F7V primary and a late M-dwarf 
secondary. The system's period is short, P = 2.230736 ± 0.000010 days, with 
an orbit indistinguishable from circular, e = 0.012 ± 0.021. We demonstrate 
generally that the surface gravity of the secondary star in a single-lined binary 
undergoing total eclipses can be derived from characteristics of the light curve 
and spectroscopic orbit. This constrains the secondary to a unique line in the 
mass-radius diagram with M/E? = constant. For HAT-TR-205-013, we assume 
the orbit has been tidally circularized, and that the primary's rotation has been 
synchronized and aligned with the orbital axis. Our observed line broadening, 
Kotsinirot = 28.9 ± 1.0 kms~^, gives a primary radius of Ra = 1-28 ± 0.04 TZq. 
Our light curve analysis leads to the radius of the secondary, = 0.167 ±0.006 
TZq, and the semimajor axis of the orbit, a = 7.54 ± 0.30 TZq = 0.0351 ± 0.0014 
AU. Our single-lined spectroscopic orbit and the semimajor axis then yield the 
individual masses, Mb = 0.124 ± 0.010 Mq and Ma = 1.04 ± 0.13 Mq. Our 
result for HAT-TR-205-013 B lies above the theoretical mass-radius models from 
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the Lyon group, consistent with results from double-lined eclipsing binaries. The 
method we describe offers the opportunity to study the very low end of the stellar 
mass-radius relation. 

Subject headings: binaries: eclipsing — binaries: spectroscopic — stars: funda- 
mental parameters — stars: low-mass — stars: rotation 



1. INTRODUCTION 



Solving for the masses and radii of stars has traditionally been accomplished through the 
analysis of double-lined eclipsing binaries, where the light of both components is detected. 
Masses and radii determined this way are fundamental and can be very accurate, because 
they rely only on Newton's laws and geometry for the analysis of the spectroscopic orbit and 
light curve, and not on models of stellar structure and evolution. In particular, analysis of the 
eclipse light curve yields the orbital inclination, and when combined with the double-lined 
spectroscopic orbit, this yields individual masses for both stars. 

There are dozens of doub le-lined ec l ipsing binaries with very accurate mass and ra- 
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NSVS01031 772 flLopez-Morales et al.M2006r) . YY Geminorum fiTorres fc RibasI l2002f ). and 
GU Bootis ( Lopez- Morales fc RibasI 2005 ). Figured] shows the observational results for these 



10 M dwarfs o n the mass-radius diagram, along with the predicted theoretical models from 
Baraffe et al.l (119981 ). All of the observed radii are larger than the theoretical predictions, 
typically by 5 to 10 percent. Another striking feature of Figure [1] is the lack of accurate 
mass-radius determinations in the lower half of the M-dwarf mass range, from CM Dra B 
(0.214 A^©), all the way down to the substellar limit near 0.075 M.q. Recently, however, the 
growing number of short-period single-lined eclipsing binaries with F- and G-star primaries 
and M-dwarf secondaries identified by wide-angle photometric surveys f or transiting planet 



promises to provi de a way to fill in this gap in the mass-radius diagram (iBouchy et al. 
Pont et all l2005al lbl. bood ) . 
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One approach to using a single-lined eclipsing binary to solve for the mass and radius of 
the unseen companion is to use stellar models together with spectroscopic and photometric 
observations of the primary to estimate its mass and radius. Then the radius ratio from the 
light curve yields the radius of the secondary, and the mass function from the spectroscopic 
orbit with the orbital inclination from the light curve can be combined to yield the mass of 
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the secondary. However, this approach is not fundamental, because it rehes on stellar models 
to characterize the primary. Tests of the mass-radius relation with the low-mass secondaries 
in such systems are no better than the validity of the models for the primaries. In particular, 
the theoretical isochrones depend on metallicity, and significant errors can result if the wrong 
metallicity is adopted. 

An alternative, and more fundamental, approach relies on the expectation that short- 
period binaries whose orbits have been circularized by tidal mechanisms must also have the 
axial rotation of both stars synchronized to the orbital period, and both rotational axes 
ali gned with the normal to the orbit, so that the two inclinations are equal, irot = "^orb (e.g. 
see lHutlll98ll : IZahrull989l ) . In this case, a measurement of the spectral line broadening due to 
rotation, V^ot sinzrot, combined with a value for Zorb from an analysis of the light curve, yields 
the actual equatorial rotational velocity, Vj-ot- If the rotation is synchronized with the orbit, 
then Prot = -Porb, and the radius of the primary can be solved. The key to this approach 
is the ability to derive accurate values for the rotational broadening of the spectral lines, 
because the primary radius can be no more accurate than the value derived for Vrot- The 
radius of the primary then sets the scale for the rest of the system, yielding the radius of the 
secondary and the semi-major axis of the orbit in actual length units, as well as the orbital 
inclination. Then Newton's revised version of Kepler's Third Law can be used to derive the 
sum of the masses, and the mass function from the single-lined spectroscopic orbit allows the 
individual masses to be determined. This approach depends on the predictions from stellar 
models only in minor ways: limb darkening coefficients are needed for the detailed analysis 
of the eclipse light curve, and the rotational broadening that is derived from the observed 
spectra can depend weakly on the metallicity that is adopted. 

The advantage of using single-lined eclipsing binaries is that it dramatically increases 
the number of low-mass stars whose masses and radii can be determined. Indeed, over the 
past 30 years only seven double-lined eclipsing systems composed of low-mass stars have been 
identified (the aforeme ntioned CM Dra, CU Cnc, N SVS01031772, YY Gem, and GU Boo, 



plus OGLE-BW3-V38 flMaceroni fc MontalbanI 120041 ) . and TrES-HerO-07621 fICreevev et al. 
20051 )). Meanwhile, more than 75 single-lined eclipsing binaries with M-dwarf secondaries 



have been discovered by wide-angle photo metric transit ing planet surveys such as Vulcan, 
TrES, and HAT in just the last five years (ILathaml 120071 ). 



Increasing the number of low-mass stars with fundamental determinations of masses and 
radii is worthwhile because of the insights these systems can yield into stellar structure. Low- 
mass stars near the hydrogen-burning limit are coo l enough that their interio r temperatures 
are on the order of the electron Fermi temperature (IChabrier fc Baraffelll997l ). causing parts 
of the stellar interior to be in the state of a partially degenerate electron gas, which means 
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a classical Maxwell-Boltzmann description of the interior does not apply. Furthermore, the 
electron number density is such that the mean inter-ionic distance is itself on the order of the 
Thomas-Fe rmi screening length, mea- ning that the electron gas is polarized by the external- 
ionic field (IChabrier &: Baraffd Il997l ). Add to this the furth er complexity that magneti c 
fields may play a role in the inner workings of low mass stars (IMullan fc MacDonaldll200ll ). 
and it becomes clear that any attempt to describe the interior of low mass stars must take 
into account the physics of both magnetic fields and partially degenerate polarized plasmas. 

Not only is it difficult to model the interiors of low-mass stars, but the usual grey 
model atmospheres are no longer applicable. The relatively low temperatures of the stellar 
atmospheres allow for the recombination of molecular hydrogen and other molecules, such 
as TiO. Therefore, accurate non-grey model atmospheres that tak e into account the effects 
of molecules must be derived and matched to the interior models (IBaraffe et al.ll2002l ). 



Thus stars near the bottom of the main sequence pose a challenge to dense matter physi- 
cists and stellar astronomers. One of the few methods of testing models for low-mass stars 
is by confronting theoretical predictions of the mass-radius relation with observations. Up 
until now there has been little data with which to constrain the possible theories. Measuring 
the masses and radii of M-dwarfs in single-lined eclipsing binaries therefore provides a new 
opportunity to test the mass-radius relation near the bottom of the main sequence. 

In this paper, we determine masses and radii for the components of HAT-TR-205- 
013, an eclipsing single -lined binary identi fied by the HAT net (Hungarian-made Automated 
Telescope Network, see lBakos et al.l (120041 )). The unseen M-dwarf secondary has a mass and 
radius of 0.124 J\4q and 0.167 TZq, with 1-a errors of 9 and 4 percent, resp ectively. This 



result places the M dwarf about 10 percent above the radius predicted by the iBaraffe et al. 



( 1l998l ) models. In this first paper we describe our analysis techniques in detail. In future 
papers we will present the results for additional M dwarf secondaries. 



2. OBSERVATIONS AND DATA REDUCTION 

2.1. HAT Photometry 

The HATnet project initiated in 2003 by G. A. B, is a wide-field survey that aims for 
the discovery of transiting planets around bright stars. It currently comprises 6 small wide- 
field automated telescopes, each of which monitors 8° x 8° of sky, typically containing 5000 



www.hatnet.hu 
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stars bright enough to permit detection of planetary transits via the typical 1% photometric 
dips they induce on their parent stars. The instruments are deployed in a two-station, 
longitude-distributed network, with four telescopes at the Fred L. Whipple Observatory in 
Arizona, and two telescopes at the Submillimeter Array at Hawaii. For a more d e tailed 



description of HAT's instrumentation, observations, and data flow, see iBakos et al.l (12002 



2004f ). 

The HAT-TR-205-013 system lies in HATnet survey field G205, centered at a = 22''55'", 
and S = -|-37°30'. 3357 observations were made of this field by the HATnet telescopes between 
5 October 2003 and 30 January 2004. Exposure times were 5 min at a cadence of 5.5 min. 
Light curves were derived by aperture photometry for the 6400 stars in G205 bright enough 
to yield photometric precision of better than 2% (reaching 0.3% in some cases). In deriving 



the l ight curves, we made use of the Trend Filtering Algorithm (IKovacs. Bakos. fc Noyes 



20051 ) to correct for spurious trends in the data. We then searc hed all the light curves for 



chara cteristic transit signals, using the Box-fitting Least Squares (IKovacs. Zucker. fc Mazeh 



2OO2I ) algorithm, which searches for box-shaped dips in the parameter space of frequency, 
transit duration, and phase of ingress. Candidate transit signals with the highest detection 
significance were then examined individually, to isolate those with the best combination of 
stellar type (preferably main-sequence stars of spectral type mid-F or later), and depth, 
shape, and duration of transit. One of these was HAT-TR-205-013, for which we identified 
a prominent periodicity with period 2.2307 days and transit depth of 0.02 mag. Figure [2] 
shows the phase-folded and fiux-normalized data for HAT-TR-205-013. 

Following our determination of HAT-TR-205-013 as a potential planetary system, we 
identified the star in the 2MASS catalog as 2MASS 23080834+3338039, which yielded J 
and K magnitudes of J = 9.691 and K = 9.408. We also found HAT-TR-205-013 in the 
Tycho catalog, as TYC 2755-36-1, with magnitudes = 11.355 and Vt = 10.729. We 
then scheduled HAT-TR-205-013 for follow-up spectroscopic observations to determine if 
the photometric dip was caused by a stellar companion. 



2.2. Follow-up Spectroscopy 



Our usual strategy for following up transiting-planet candidates identified by wide-field 
photometric surveys is to start with an initial spectroscopic reconnaissance, to see if there 
is evidence for a stellar companion th at is responsib le for the observed light curve. We have 
used the CfA Digital Speedometers (ILathaml Il992l ) on the 1.5-m Wyeth Reflector at the 
Oak Ridge Observatory in the Town of Harvard, Massachusetts and on the 1.5-m Tillinghast 
Reflector at the Fred L. Whipple Observatory on Mount Hopkins, Arizona to obtain single- 



- 6 - 



order echelle spectra in a wavelength window of 45 A centered at 5187 A with a resolution of 
8.5 kms~^ and a typical signal-to-noise ratio per resolution element of 15 to 20. For slowly 
rotating solar-type stars these spectra deliver radial velocities accurate to about 0.5 kms"^, 
which is sufficient to detect orbital motion due to companions with masses down to about 5 
or 10 J for orbital periods of a few days. 

Spectroscopic follow-up observations of transiting-planet candidates are easier to sched- 
ule than photometric observations, because the radial velocity varies throughout the entire 
orbital phase, while the transit light curve has a duty cycle of only a few percent, e.g. 3 
hours out of 3 days. Spectroscopy has the additional advantage over photometry that we 
can use our spectra to classify the parent star, deriving effective temperature, rotational 
velocity, and surface gravity by correlating the observed spectra with a library of synthetic 
spectra. This information is often useful in rejecting host stars that are too hot, or are 
rotating too rapidly, or are evolved (in which case we assume the giant must be the bright 
star in a blended triple). Thus at CfA we usually start with a spectroscopic reconnaissance 
before attempting follo w-up photometry or very precise radial-velocity observations (e.g. see 
O 'Donovan et ahlboOTl ). 



Starting in 1999 with transiting-planet candidates provided by the Vulcan team (IBorucki et al. 



200ll ). we quickly learned that the vast majority o f the planet can didates were actually eclips- 
ing binaries masquerading as transiting planets (jLathanj|2003l ). One of the most common 
imposters was an F- or G-star primary eclipsed by a late M-dwarf secondary. Such systems 
produce light curves similar to planets transiting solar-type dwarfs, because M dwarfs near 
the bottom of the main sequence have roughly the same size as giant planets. Nevertheless, 
they are easy to distinguish from planets, because their masses are two orders of magnitude 
larger, and the reflex orbital motion they induce in their parent stars has an amplitude of 
at least several kms~^. 

We determined the rotational and radial velocities of HAT-TR-205-013 by cross corre- 
lation of the observed spectra against templates dra wn from a library of synthetic spectra 
calculated by Jon Morse for a grid of iKuruca (119921 ) stellar atmospheres. The library grid 
has a spacing of 250 K in effective temperature, T^s] 0.5 in log surface gravity, logg; and 
0.5 in log metallicity relative to the sun, [Fe/H]. For the final radial-velocity determinations 
we adopted the template with the Tes, ^ogg, and Vmt values that gave the highest value for 
the peak of the correlation coefficient, averaged over all the observed spectra, assuming solar 
metallicity. For the correlation analysis we used xcsao (IKurtz fc Minklll998l ) running inside 
the IRAFli environment. 



^IRAF (Image Reduction and Analysis Facility) is distributed by the National Optical Astronomy Ob- 
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Our very first spectrum of HAT-TR-205-013 revealed that the spectral lines of the 
primary were broadened by about 30 kms~^ of rotation. In our experience this is a strong 
suggestion that the companion is a star with enough mass to synchronize the rotation of 
the primary with the orbital period. The second exposure showed that the velocity varied 
by several kms~^. Additional spectra soon revealed a spectroscopic orbit with a period 
that matched the photometric period and an orbital semi-amplitude, K, of about 20 kms^^ 
implying that the secondary was a small M dwarf. Altogether we accumulated 23 spectra 
yielding the radial velocities reported in Table [H which were calculated using a template 
with Teff = 6250 K, logs' = 4.0, [Fe/H]=0.0, and Vi-ot = 30 kms~^. The parameters of our 
orbital solution for HAT-TR-205-013 are reported in Table [2] and the corresponding velocity 
curve and observed velocities are plotted in Figure [31 

Because the value of the rotational velocity is critical for our determination of the radius 
of the primary star, we evaluated both the internal precision of our determinations and 
possible systematic errors resulting from uncertainties in the surface gravity and metallicity 
of the primary star. Our final results for the mass and radius of the primary star imply a 
surface gravity logs' = 4.24, halfway between the nearest templates in our library of synthetic 
spectra. Therefore we evaluated the rotational velocity at both logs = 4.0 and 4.5, and also 
at two metallicities, [Fe/H] = 0.0 and —0.5. At each of the four combinations of logs s-iid 
[Fe/H] we ran correlations over a wide range of temperatures, 5500 to 7250 K, and rotational 
velocities, 10 to 50 kms~^. At each gravity and metallicity we interpolated to find the 
temperature and rotation that gave the highest value of the correlation coefficient averaged 
over all the observed spectra weighted by the photon statistics. For these experiments we 
did not use two of our observed spectra that had much lower exposure levels than the other 
21 observations. The results are reported in Table [31 They show that the accuracy of our 
rotational- velocity determination at a given logs ^"^^ [F^/H] is not seriously limited by the 
scatter from the individual exposures, despite the relatively low signal-to-noise ratio of our 
observed spectra. The uncertainty in the mean rotational velocity at each logs ^"^^ [F^/H] 
due to internal scatter is less than 0.3 kms~^. 

The results reported in Table [31 also show that the systematic errors due to uncertainties 
in the gravity and metallicity are not serious. For example, the dependence of rotational 
velocity on metallicity is rather weak, only 0.2 kms~^ for a change of 0.5 in [Fe/H]. On the 
other hand, the dependence on gravity is rather strong, about 1.0 kms~^ for a change of 0.5 
in logs- Fortunately, the uncertainty in the actual gravity of the primary star is quite small, 
less than 0.2 in logs, so the rotational velocity interpolated between the two gravities should 



servatories, which are operated by the Association of Universities for Research in Astronomy, Inc., under 
contract with the National Science Foundation. 
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be good to perhaps 0.2 kms~^. Unfortunately we do not yet have a direct measure of the 
metalhcity of HAT-TR-205-013, but the fact that the templates with [Fe/H] = 0.0 give a 
better match to the observed spectra than the templates with [Fe/H] = —0.5 suggests that the 
metallicity is probably within 0.5 dex of solar. In the absence of an accur ate metallicity we 



adop t a value of [Fe/H] = —0.2, which is typical of the solar neighborhood (INordstrom et al. 



20041 ). The temperature corresponding to this metallicity that we derive from our spectra is 
Teff = 6295 K, but with a systematic error that could exceed 200 K because of the degeneracy 
between temperature and metallicity in the analysis of our spectra. The alternative approach 
of deriving a temperature from photometric indices is hampered by the lack of accurate 
photometry and the possibility that there may be significant reddening. The distance that 
we derive from the Stefan-Boltzmann law using R = 1.27 ± 0.04 TZq, Tcs = 6295 ± 200 K, 
Vt = 10.7 29 =b 0.048 (w hich corresponds to ^ = 10.67 ± 0.05), and bolometric correction 



—0.011 (lFlowerlll996l ) is 232 ± 18 pc. From our sp ectroscopi c temperature and logg, we 



estimate the spectral type of the primary to be F7V (ICoxll2000l ). 



Thus our experiments with different gravities and metallicities suggest that the system- 
atic errors due to uncertainties in the template parameters may be less than 0.4 kms~^. 
When combined in quadrature with the internal precision estimate of 0.3 kms~^, this sug- 
gests that the total error in our final interpolated rotational velocity could be as small as 0.5 
kms~^. Previous experience however, would indicate that the actual uncertainty is higher 
than this internal estimate. Therefore, to be conservative we adopt an uncertainty of 1.0 
kms~^ and use l^otsinvt = 28.9 ± 1.0 kms~^ for our determination of the radius of the 
primary and subsequent analysis. 

One could also determine the value of Kotsinirot by sp ectroscopically measu ring the 



amplitude of the Rossiter-McLaughlin Effect during echpse (IGaudi fc Wind 120071 ) . In the 
case of HAT-TR-205-013 the amplitude of this effect would be about 0.5 km s^^, too small to 
be measured with the CfA Digital Speedometers, but potentially observable with an accuracy 
of a few percent by instruments such as HIRES on Keck 1. 

Our approach relies on the assumption that the stellar rotation has been synchronized to 
the orbital period. One way to test this assumption might be to use very precise photometry 
to derive a rotational period for the star. We did look for sinusoidal variations in the HAT 
photometry but did not find anything significant near the orbital period. We did detect a 
marginally significant variation with half the orbital period and amplitude of about 1 mmag, 
which is consistent with the the expected ellipsoidal deformation of the primary star. A 
second, more global test would be to use stellar models to estimate the mass and radius of 
the primary, but a prerequisite for such an analysis would be accurate determinations of the 
temperature and metallicity, presumably from high-quality spectra, and such results are not 
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yet available. 



2.3. Follow-up KeplerCam Photometry 

To provide a high-quality light curve for the analysis of the primary eclipse of HAT- 
TR-205-013 we used KeplerCam on the 1.2-m telescope at the Fred L. Whipple Observatory 
on Mount Hopkins, Arizona. KeplerCam utilizes a monolithic 4Kx4K Fairchild 486 CCD 
that provides a 23' field and a pixel size of 0.34". We used the predicted eclipse times from 
our spectroscopic orbit to schedule observations on the night of 22-23 October 2005. We 
successfully observed a full eclipse, alternating between the Sloan g and i bands. 

The seeing was approximately 2" FWHM throughout the night, but we deliberately 
defocused the telescope to get images with 3" FWHM in order to keep the peak counts 
in the images well below saturation. Over the course of the observations the focus of the 
telescope, which depends weakly on temperature, was adjusted three times to keep the image 
size near 3" FWHM. Because we used automatic guiding, the centroid of the images moved 
less than 3 pixels over the duration of the observations. For readout we binned the pixels 
2x2, which gave a total readout time including overhead of 12 seconds. Exposure times 
were 30 seconds for the g band and 10 seconds for the i band. The images were obtained in 
sequences of 3 g-hand exposures followed by 6 i-band exposures. All told, we were able to 
collect a total of 297 images in the g band and 588 in the i band. Some thin cirrus clouds 
were present at the beginning of the night, with no noticeable degradation of the light curves. 
A quarter Moon rose during the observations after the end of egress, and this contributed 
to a slight increase in photometric scatter after egress. 

All the images were reduced by applying an overscan correction and then subtracting 
the two-dimensional residual bias pattern. After correcting for shutter effects, we flattened 
each image using a normalized set of combined twilight images. To produce the light curve, 
we used the first image from each filter in an observing sequence as our astrometric reference 
for identifying the same stars in subsequent images. We then determined the relative shift 
between images to relocate each star in the following images. We measured the flux of each 
star in a 6.7" circular apertur e around the po sition derived from the astrometric fit using 



daophot/phot within IRAF (jTodylll986l . Il993l ) . We estimated the sky in an annulus around 



each star with inner and outer radii of 9.4" and 13.4" using the sigma-rejection mode. 

We iteratively selected comparison stars by removing any that showed unusual noise or 
variability in their differential light curves, obtaining 21 comparison stars for the g band and 
37 for the i band. Based on a weighted mean flux from the comparison stars, we calculated 
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an extinction correction, and then applied this to each comparison star. The typical rms 
residual for both bands was 1.5 mmag, a value that increased slightly toward the end of the 
night due to the increased sky brightness from the rising Moon. The main contributor of noise 
to the photometry from individual images was atmospheric scintillation, which accounted 
for over 65% of our calculated rms (jYounpl967l ). The individual photometric measurements 
in the g and i bands are reported in Tables H] and [5l respectively. 



3. LIGHT-CURVE ANALYSIS 

Our KeplerCam light curves for HAT-TR-205-013 give good coverage of the eclipse 
centered at HJD 2453666.747 ± 0.001, with more than two hours of coverage both before the 
start of ingress and after the end of egress. The eclipse itself lasts about 3 hours between 
first and fourth contact and is about 2 percent deep, in agreement with the discovery light 
curve from HAT. Our KeplerCam light curves clearly show the effects of limb darkening on 
eclipse shape as a function of wavelength: the i-band light curve is slightly shallower and 
possesses a flatter bottom than the g-hand curve. The portions of the light curves before 
and after eclipse showed a slight drift, which we removed with a linear fit. The resulting 
light curves are plotted in Figure HI 

From our fits to the KeplerCam light curves we derive transit centers of 2453666.7465 ± 
0.0005 and 2453666.7473 ± 0.0005 in the g and i bands, respectively. When we com- 
bine these transit times with the results from the HAT photometry, we get the ephemeris 
2453666.74748 ± 0.00018 + A x 2.230736 ± 0.000010. 

Knowing the duration and depth of the eclipse, together with the orbital period, we 
were able to derive rough values for the ratio of semimajor axis and radius of the secondary 
to the radius of the primary star using the relations 



a I f P 



Ra n- V Attr 



'1^ 



^^VS (2) 

Ka 

where a is the orbital semimajor axis, -Ra and Rb are the primary and secondary stars' radii, 
P is the orbital period, Attr is the transit length, and 6 is the transit depth in relative flux. 
For our values of P ~ 2.23 days, Attr ~ 2.5 hours and 6 ~ 0.02, we obtained a/RA ~ 5.93 
and Rb/Ra ~ 0.133 from this first approximation. 

We then built a two dimensional grid of light-curve fits in q/Ra and R-q/Rai centered 
on the rough values obtained above. The impact parameter h was introduced as a third 
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dimension in the grid, varying from 6 = (central transit) to 6 = 1 + Rb/Ra (grazing 
transit). The grid steps were 0.002 in g/Ra, 0.00001 in R^/Ra, and 0.001 in b. We generated 
synthetic hght-curves for each combination of parameters using the routines provided by 



Mandel fc Agoll (|2002| ). together with the quadratic hmb-darkening coefficients derived by 
Claretl (2004), using the temperature and surface gravity for the primary we had derived 
spectroscopically. We adopted solar values for the metallicity and surface turbulence of the 
primary star. The exact coefficients we used were Ui = 0.4238 and U2 = 0.3250 in the g, 
and Ml = 0.1814 and U2 = 0.3723 in the i band. Uncertainties in our values of the limb- 
darkening coefficients, largely a result of the uncertainty in the spectroscopic temperature 
measurement, had a negligible effect on our final results. 

To identify the synthetic light curve that gave the best fit to the actual data, we looked 
for the fit with the smallest value of x^? inspecting the contours to ensure that we had 
found the global minimum. We also inspected plots of the fits to the data as a qualitative 
check. We applied this procedure to both the g and i light curves. The results are summarized 
in Table [61 

To calculate the uncertainties for our values of a/R^, Rb/Ra, and b, we performed 
additional fits wherein we fixed one of the parameters at a value differing slightly from the 
best-fit value, and fit for the other two parameters. We changed this fixed parameter until 
the eventual fit achieved a that was 1-a away from our best fit parameters' x^- We 
then used the difference between this 1-a value and our best fit value as the uncertainty 
for that particular parameter. We examined the correlation between the three different 
parameters by making a final fit to the data using the Levenberg -Marquardt meth od - a 



combination of the Inverse-Hessian and Steepest Descent methods (jPress et al.lll992l ). This 
allowed us to compute the covariance matrix for the three fit parameters, and thus the 
correlation coefficients for the parameters. For all three, the correlation with any of the 
other two was always smaller than 0.5 in both bands. While this is not negligible, we found 
that the contribution to the total estimated error was minor. Of course, the Levenberg- 
Marquardt method can also be used to calculate the uncertainties in a/i?A, Rb/Ra, and b. 
However, the uncertainties that we derived from our grid analysis were substantially larger. 
To be conservative, we adopted the grid uncertainties instead of the Levenberg-Marquardt 
uncertainties. Figure shows contour plots of for our fits. Table E] reports the final fitting 
results and errors, and Table [7] lists the correlation coefficients. The KeplerCam data and 
best-fit light curves are shown in Figure |H Simultaneous fitting of the light curves in both 
bands produced fits with inferior x"^ statistics, a result of the slight difference in Rb/Ra in 
the two bands. We therefore use a simple average of the light curve results as our adopted 
values for the parameters in Table [61 
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3.1. MASSES AND RADII FOR HAT-TR-205-013 

Using the values that we have measured from the echpse hght curves, as well as the 
observed spectroscopic orbit parameters, we may restrict the location of HAT-TR-205-013 
A and B on the mass-radius diagram to unique curves described completely by these observ- 
ables. To emphasize this, we have written the observable values in brackets in the following 
derivation. 

To begin, we use Newton's revised version of Kepler's Third Law and the spectroscopic 
mass function from the orbital solution: 

where [Ka] is the semi-amplitude of the spectroscopic orbit for the primary in km s^^, [e] 
the eccentricity of the orbit, and G is the Gravitational Constant. With these two equations, 
the two unknowns (Ma and Mb) may be found: 



" G [PP y 2x<jsmi„,t I 

M, = ^ J^V^W? (6) 

G [P] smzorb 

Both the orbital inclination and the semimajor axis can be expressed in terms of the 
observables [a/i?A], [Rb/Ra] and [b], and the unknown stellar radii Ra and Rb- 

By taking our value for the impact parameter [b], which measures the projected separa- 
tion of the primary and secondary stars at the mid-point of the eclipse, the orbital inclination 
can be expressed in terms of the observables [a/i?A] and [b]: 

[6] = [a/i?A]cosiorb (7) 

/ [612 \l/2 

= V-Wr^^) 

Similarly, the semimajor axis may be written in terms of the observables [a/ Ra] , [-Rb / Ra] 
and [b], and the unknown stellar radii Ra and Rb' 



(9) 
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Substituting these values into Eq.(5) and (6) gives expressions for the mass of each 
component as functions of only the observables and the respective stellar radii: 



Ma - ^ [a/R.r 1 '^^-J ^\-^'^' I Rl (10) 

G[P? V 27r(l - [b]y[a/RA?f' [a/i?A] i?A ' 



27r / [a/RA] V [K^] y^l - [e] 
G[P] [[Rb/Ra]J il-[b]y[a/RA?f' 



Mb = ro'7o\ ) -^B (11) 



Therefore, solely from what we are able to measure using the eclipse light curves and 
spectroscopy, we may confine HAT-TR-205-013 B to a single curve on the mass-radius dia- 
gram that goes as Mb/-R| = constant. This constant is di rectly proportiona l to th e surface 



gravity of the object, which has been noted previously by ISouthworth et al.l (120041 ) . In our 
case, the quality of our photometry allows us to measure this constant extremely accurately; 
the uncertainty region around this surface gravity line is on the order of the line width itself 
when it is plotted. 

Through the assumed synchronization of the primary's rotation to the orbital period we 
may locate HAT-TR-205-013 B on this gravity determination curve, but it is important to 
note that the line itself is defined without having made any assumptions about the system. 
Indeed, it is possible to calculate a similar gravity curve for any eclipsing system, or for a 
system containing a transiting planet. All that is required is a good-quality light curve and 
spectroscop i c orbi t. This has previously been done for the case of a transiting planet by 



Winn et all (120061 ). 



To place HAT-TR-205-013 on the mass-radius diagram more specifically, we note that 
that the measured eccentricity of our spectroscopic orbit for HAT-TR-205-013 is indistin- 
guishable from circular, and therefore, for the reasons described in the Introduction, we 
assume that the spin axes of both stars have been aligned with the orbital normal and that 
the rotation of both stars has been synchronized to the orbital period. This allows us to use 
the observed rotational line broadening of the primary to solve for the radius of the primary 
in linear units, which in turn allows us to convert the orbital size and secondary radius into 
linear units from the values of [a/Rp^] and [Rb/Ra] derived from the light curves. 

Using the assumption of synchronization, and that iorb = "^rot, we see by inspection that 

R^ = (12) 

27r sin i^rh 

R. = g "^"'^'"'-' ifie/iiAl (13) 

/vr sm torh 

where [VJ-ot sin ijot] is the projected rotational broadening of the primary derived from its 
observed spectra. We may now substitute in Eq.(8) for siniorb to get both radii in terms of 
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our observables: 



Ra = — TTT [Kotsinirot] (14) 

2n{l-[b]y[a/RA]'f' 



27r(l-[6]V[a/i?A] 



= ,J,,!^ ,r. 12^1/2 [-^B/i?A] [VrotSin^,ot] (15) 



By combining these two statements with Eq.(9) and (10), we arrive at expressions for 
the masses of each component in terms of just the observable quantities: 

27rG(l_[6]2/[a/i?^]2)3/2 \^ [a/i?A][KotsinZrot]; ^ ^ ^ ^ 

= ^ ,^ ,!!(;^^^^ \Kk\ v^W^ [Kot sin^^ot]^ (17) 



27rG (1 _ [6]2/[a/i?A 



The results for the masses and radii for both components of HAT-TR-205-013 are pre- 
sented in Table [H The errors were estimated using Monte-Carlo simulations and were com- 
pared with the results of formal error propagation, including the correlation coefficients 
derived from the light-curve fits. Both approaches delivered similar results. The mass and 
radius obtained for the primary star are essentially the same for both the g and % light curves, 
but the mass and radius for the secondary differ by 0.8 and 3 percent, respectively. This 
radius difference between the two light curves is close to l-cr, and may be due to uncertainties 
in the limb-darkening coefficients. Our adopted values are based on the average values of 
the light curve parameters. 



4. DISCUSSION 



In Figure [6] we plot our mass and radius for the M-dwarf secondary HAT-TR-205- 
013 B on a rnass-ra dius diagram, together with isochrones for ages of 0.5 and 5 Gyr from 
Baraffe et al.l (119981 ). We also plot the res ults for 11 M dwarf s econdaries f r om the s ample 
of OGLE planetary candidates analyzed by Bouchy et aL ( 2005 ): Pont et al. ( 2005a| jbl. I2OO6 ) 
and listed in Tabl e El For the systems OGLE-TR-34 JBoucIiv et al.ll2005h. OG LE-TR-120 
dPont et al.ll2005bh. andthe low mass systems OGLE-TR-122 JPont et al.ll2005ah and OGLE- 
TR-123 (jPont et al.ll2006l ) the authors had to use stellar models to to estimate the masses 
and radii of the primaries without the assumption of synchronization, as synchronization 
implied masses and radii that were inconsistent with the spectroscopic observations. For the 
other seven systems, they were able to assume synchronization and to derive the radius of the 
primary from the observed rotational line broadening. In general the agreement between the 
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OGLE results and the lBaraffe et al.l (Il998l ) isochrones looks promising, but the observational 
uncertainties are still too large to allow a critical test of the theoretical models. The OGLE 
systems are all much fainter than HAT-TR-205-013, which presents significant challenges 
for both the spectroscopic and photometric follow-up observations. Spectroscopy with the 
resolution and signal-to-noise ratio suitable for determining accurate values for rotational 
broadening requires time on large telescopes, and photometry for high-quality light curves 
also requires large telescopes to achieve the needed photon statistics. Eclipsing binaries 
identified by wide-angle surveys are much brighter and therefore less challenging on both 
counts. 

Our value for the radius of the M-dwarf secondary in HAT-TR-205-013 lies 11 percent, 
about 3-0", above the theoretical isochrones. This divergence is further reinforced by Eq.(ll), 
which, as has been previously noted, restricts the position of HAT-TR-205-013 B to lie on a 
single line that is determined by the surface gravity of the object. This gravity curve does 
not rely on any prior assumptions about the HAT-TR-205-013 system, nor does it depend 
upon our measured value of V^t sinirot, which is the biggest contributor of uncertainty to our 
final results. We use the assumption of synchronization and the spectroscopically measured 
V^rotsinirot to place HAT-TR-205-013 B at a specific location along the curve, but it is 
important to note that in the region that we find HAT-TR-205-013 B, the curve of allowable 
locations runs nearly parallel to the theoretical models. This is illustrated in Figure [6] by 
the red line that passes through our point for HAT-TR-205-013 B. Thus the conclusion that 
the theoretical models predict a radius for HAT-TR-205-013 B that is too small by about 10 
percent is on much firmer ground than the error bar in the observed radius might suggest. 
Indeed, it would require a 6-a difference in V^rot sinirot to place HAT-TR-205-013 B onto the 
Baraffe models. 

Our result for HAT-TR-205-013 B supports the suggestion from the results for double- 
lined eclipsing binaries plotted in Figure [T] that the models predict radii for M dwarfs that 
are too small by up to 1 percent. This discrepan cy has been n oted before, for example 



by [Torres &: RibasI (120021 ) in the case of YY Gem. [Torres et al.l (120061 ) raised the issue of 
whether short-period eclipsing binaries are representative of isolated field stars and wide 
binaries where tidal forces are negligible. They suggested that the rapid rotation of the 
stars in these systems caused by tidal synchronization might give rise to enhanced magnetic 
activity, thus decreasing the efficiency of energy transport in the convective envelopes and 
leading to infla t ed ste llar radii. For low mass stars, this effect is examined in more detail by 



Lopez-Moraled (120071 ). 



In the case of HAT-TR-205-013, we see no evidence in the photometry of star spots on 
the primary star, which would be tell-tale indicators of enhanced stellar magnetic activity. 
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Though HAT-TR-205-013 A is rapidly rotating, the lack of magnetic ativity is not suprising, 
given its spectral type (F7). The star's outer convective layer is relatively shallow, and 
i t is not unusual fo r rapidly rotating stars of this type to lack strong magnetic activity 



flTorres et al.ll2006f ) 



In some instances it may be possible to independently determine the rotational period of 
the primary through high-quality light curves used to definitively identify photometric varia- 
tion outside of eclipse. This would serve as a check to the assumption of tidal synchronization 
in the system. 

In future papers we will present the results for several additional single-lined eclipsing 
binaries with circularized orbits. 
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Table 1. Individual Radial Velocities 



HJD 








(days) 


(kms~ 




(kms~^) 


2453034.45642 


-2. 


,02 


1.38 


2453035.47574 


-10, 


.11 


1.61 


2453035.58018 


-18, 


.52 


1.01 


2453036.48778 


-12, 


.93 


1.53 


2453037.46565 


-0. 


,47 


1.43 


2453037.61215 


-6. 


,83 


0.91 


2453038.46413 


-25, 


.72 


1.48 


2453038.57874 


-19, 


.76 


1.15 


2453040.47360 


-28, 


.58 


1.24 


2453042.58686 


-27, 


.79 


0.91 


2453043.58338 


+4, 


.84 


1.23 


2453044.58422 


-19 


.73 


1.14 


2453045.57911 


-6, 


,42 


0.73 


2453046.46373 


-2. 


,11 


0.80 


2453046.60000 


-10, 


.47 


0.77 


2453047.50881 


-20, 


.45 


1.49 


2453047.58731 


-14, 


.83 


0.95 


2453543.94910 


-4. 


,01 


1.10 


2453658.69572 


-20, 


.53 


1.16 


2453659.75967 


+3, 


.14 


2.28 


2453659.78398 


+2, 


.85 


1.09 


2453660.70213 


-25, 


.98 


1.57 


2453664.70202 


-19, 


.39 


1.21 
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Table 2. Spectroscopic Orbital Parameters 



Parameter 


Value 


P (days) 


2.23072 ± 0.00005 


7 (kms~"^) 


-9.83 ±0.30 


K (kms-^) 


18.33 ±0.47 


e 


0.012 ±0.021 




143 ± 90 


Epoch (HJD) 


2, 453, 198.61 ±0.56 


Nobs 


23 


— C rms (kms~^) 


1.06 


/(M) {Me') 


0.00142 ± 0.00023 


a A sin i (Gm) 


0.562 ± 0.030 
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Table 3. Rotational Velocity Results 



log5,[Fe/H] 




< Vrot > 


(t{< Vrot >) 


Correlation 




(K) 


(kms~^) 


(kms~^) 


Coefficient 


4.0,0.0 


6340 


29.4 


0.25 


0.826 


4.5,0.0 


6540 


28.4 


0.24 


0.823 


4.0,-0.5 


5960 


29.2 


0.21 


0.821 


4.5,-0.5 


6150 


28.2 


0.24 


0.816 


Adopted: 










4.24,-0.2 


6295 


28.9 


1.0 





-22- 



Table 4. g Band Photometry 



HJD Flux 



2453666.575985 1.00054 

2453666.576472 0.99856 

2453666.576946 1.00108 

2453666.579226 1.00084 

2453666.579712 1.00008 

2453666.580198 0.99985 

2453666.582501 0.99998 

2453666.582976 0.99882 

2453666.583474 1.00159 



Note. — Table 4 is pre- 
sented in its entirety in the 
electronic edition of the As- 
trophysical Journal. A por- 
tion is shown here for guid- 
ance regarding its form and 
content. 

column (1): Heliocentric 
Julian Date, 

column (2): Normalized in- 
strumental flux. 
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Table 5. i Band Photometry 



HJD Flux 



2453666. 


,574226 


0, 


.99951 


2453666, 


,574469 


0, 


.99803 


2453666. 


,574724 


1, 


.00096 


2453666, 


,574967 


0, 


.99846 


2453666, 


,575233 


0, 


.99628 


2453666, 


,575488 


1, 


.00014 


2453666, 




1 


.00062 


2453666, 


,577698 


0, 


.99886 


2453666, 


,577965 


0, 


.99941 


2453666, 


,578208 


0, 


.99613 


2453666, 


,578474 


0, 


.99957 


2453666, 


,578728 


0, 


.99962 


2453666, 


,580719 


0, 


.99837 


2453666, 


,580974 


0, 


.99976 


2453666, 


,581228 


0, 


.99761 


2453666, 


,581494 


1 


.00025 


2453666, 


,581772 


1 


.00113 


2453666, 


,582027 


0, 


.99823 



Note. — Table 5 is pre- 
sented in its entirety in the 
electronic edition of the As- 
trophysical Journal. A por- 
tion is shown here for guid- 
ance regarding its form and 
content. 

column (1): Heliocentric 
Julian Date, 



-24- 



column (2): Normalized in- 
strumental flux. 
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Table 6. Light-Curve Fit Results 



Parameter 


g Band 


i Band 


Adopted 


a/RA 


5.93 ±0.15 


5.91 ±0.16 


5.92 ±0.11 


Rb/ Ra 


0.1330 ±0.0010 


0.1288 ±0.0007 


0.1309 ±0.0006 


b 


0.36 ±0.06 


0.37 ±0.07 


0.365 ±0.046 
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Table 7. Correlation Coefficients 



Coefficient g Band i Band 



{a/RA,RB/RA) 0.28 0.27 
ia/RA,b) -0.21 -0.42 
{RB/RA,b) 0.04 0.01 
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Table 8. Physical Parameters for HAT-TR-205-013 



Parameter 


g Band 


i Band 


Adopted 


Ra {Uq) 


1.04 ±0.14 
1.28 ±0.04 


1.03 ±0.14 
1.28 ±0.04 


1.04 ±0.13 
1.28 ±0.04 


Mb [Mq) 

Rb (7^0) 


0.124 ±0.011 
0.169 ±0.006 


0.123 ±0.011 
0.164 ±0.006 


0.124 ±0.010 
0.167 ±0.006 


a (AU) 


0.0351 ± 0.0015 


0.0351 ± 0.0015 


0.0351 ± 0.0014 
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Table 9. MASSES AND RADII FOR Low-Mass Stars 



Name M (Mq) R {Uq) Type Ref. 



OGLE-TR-123 B 


0.085 


± 





oil 


0.133 


± 





009 


SBl 


EB 


1 


OGLE-TR-122 B 


0.092 


± 





009 


0.120 


± 





018 


SBl 


EB 


2,3 


OGLE-TR-106 B 


0.116 


± 





021 


0.181 


± 





013 


SBl 


EB 


3 


HAT-TR-205-013 B 


0.123 


± 





oil 


0.167 


± 





007 


SBl 


EB 


13 


OGLE-TR-125 B 


0.209 


± 





033 


0.211 


± 





027 


SBl 


EB 


3 


CM Dra B 


0.2136 


± 





0010 


0.2347 


± 





0019 


SB2 


EB 


4,5 


CM Dra A 


0.2307 


± 





0010 


0.2516 


± 





0020 


SB2 


EB 


4,5 


OGLE-TR-78 B 


0.243 


± 





015 


0.240 


± 





013 


SBl 


EB 


3 


OGLE-TR-5 B 


0.271 


± 





035 


0.263 


± 





012 


SBl 


EB 


6 


OGLE-TR-7 B 


0.281 


± 





029 


0.282 


± 





013 


SBl 


EB 


6 


OGLE-TR-6 B 


0.359 


± 





025 


0.393 


± 





018 


SBl 


EB 


6 


OGLE-TR-18 B 


0.387 


± 





049 


0.390 


± 





040 


SBl 


EB 


6 


CU Cnc B 


0.3890 


± 





0014 


0.3908 


± 





0094 


SB2 


EB 


7 


OGLE-BW3-V38 B 


0.41 


± 





09 


0.44 


± 





06 


SB2 


EB 


8 


CU Cnc A 


0.4333 


± 





0017 


0.4317 


± 





0052 


SB2 


EB 


7 


OGLE-BW3-V38 A 


0.44 


± 





07 


0.51 


± 





04 


SB2 


EB 


8 


OGLE-TR-120 B 


0.47 


± 





04 


0.42 


± 





02 


SBl 


EB 


3 


TrES-HerO-07621 B 


0.489 


± 





003 


0.452 


± 





050 


SB2 


EB 


9 


TrES-HerO-07621 A 


0.493 


± 





003 


0.453 


± 





060 


SB2 


EB 


9 


NSVS01031772 B 


0.4982 


± 





0025 


0.5088 


± 





0030 


SB2 


EB 


10 


OGLE-TR-34 B 


0.509 


± 





038 


0.435 


± 





033 


SBl 


EB 


6 


NSVS01031772 A 


0.5428 


± 





0027 


0.5260 


± 





0028 


SB2 


EB 


10 


YY Gem A & B 


0.5992 


± 





0047 


0.6191 


± 





0057 


SB2 


EB 


11 


GU Boo B 


0.599 


± 





006 


0.620 


± 





020 


SB2 


EB 
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Fig. 1. — The mass-radius diagram for 10 stars in 5 double- lined eclipsing binaries each 
composed of two M dwarfs, and with errors better than 3 percent. 
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2.— The phase-folded HATnet hght curve for HAT-TR-205-013. 
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Fig. 3. — The velocity curve for our orbital solution for HAT-TR-205-013, together with the 
individual observed velocities. The lower panel shows the 0-C velocity residuals from the 
orbital solution. 
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Fig. 4. — KeplerCam light curves for HAT-TR-205-013 in the SDSS g and i bands. Contin- 
uous hnes show the best fit synthetic hght curves for each. 
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Fig. 5. — Contours of for the results from fits to the hght curves in the 51-band (left 
panels) and i-band (right panels). For each band the three panels show the projections onto 
the three possible planes involving b, a/R\, and Rb/Ra- The l-cr, 2-a, and 3-a contours are 
plotted. 
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Fig. 6. — T he mass-rad i us diagram for M dwarfs in single-lined eclipsing binaries. The M 
dwarfs from iPont et al.l (l2005al Jbl. l2006l ) are plotted as open circles. The red line passing 
through the point for HAT-TR-205-013 B shows the constraint imposed on its location by 
Eq.(ll) and our observed quantities, without making any explicit assumptions (such as 
synchronization) about the system. Assuming synchronization, the hash marks on the line 
show the effect that differences of ± 1, 2, & 3 kms~^ in V^ot have on our final results. 



